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GÉNÉRALISATION DES THÉORIES DE LA CONVECTION CALORIFIQUE. 
NOTION DE CONVECTION VIVE, CONSÉQUENCES TECHNIQUES. I 
Par Yves ROCARD et MARCEL VÉRON. 
Sommaire. — Le présent travail a d’abord pour objet d'approfondir et d’étendre la notion de convection 


calorifique, en vue de permettre un calcul approché des coefficients de convection dans des cas moins 
spéciaux que ceux actuellement pris en considération, puis d’inspirer une meilleure appréciation, 
peut-être aussi — dans certains cas — une conception plus rationnelle des matériels thermiques. 
I] précise, commente et généralise trois communications récemment présentées à l’Académie (1). 

Nous serons amenés à distinguer deux sortes de convections : la convection morte, forme classique 
où le fluide chauffant ou refroidissant ést thermodynamiquement passif, et la convection vive, 
nouvelle forme qui apparaît dans les cas où le fluide participe aussi à la création ou à la disparition de 
la chaleur sensible qu’il véhicule; tel est par exemple le cas d’une flamme, mélange gazeux en cours de 
combustion. Les réactions, la diffusion, le rayonnement, et toutes leurs combinaisons possibles, 
introduisent autant de variantes de la convection vive. 

Pour la convection accompagnée d’une réaction, on présentera le coefficient de convection comme 
formé du coefficient de convection habituel, plus un ferme correctif que nous appellerons terme vif; 
celui-ci devient sans objet quand le fluide reste passif. 

On se tromperait très grossièrement si l’on appliquait le coefficient de convection morte à l’esti- 
mation du pouvoir chauffant d’une flamme venant lécher la sole d’un four, le terme correctif 
attaché à la convection vive étant alors non seulement très important, mais souvent largement 
prépondérant, au moins pour les flammes dites « actives ». 

_Il fait jouer des facteurs tout à fait étrangers à la convection classique mais familiers aux techniciens 
intéressés, telle la quantité de chaleur {qui se trouve créée ou détruite dans l’unité de volume, 
par unité de temps; aussi sa considération éclaire-t-elle d’une lumière nouvelle une foule de faits 
plus ou moins connus, faits qu’on s'était peu soucié d’expliquer jusqu'ici. 

Ainsi apparaît tout l'intérêt pratique des nouvelles notions. 

Pour les rendre strictement quantitatives dans les différents cas d’espèce, beaucoup d’eflorts 
conjugués seront nécessaires; mais elles semblent dès maintenant de nature à réaliser l’harmonie 
entre les données empiriques, les développements théoriques et le bon sens — à leur commun bénéfice. 


I. — L’équation générale de la chaleur 
et le coefficient de convection. 


1. L'équation restreinte de la chaleur. 
L’équation classique qui règle la distribution de la 
température dans un fluide exprime que la quantité 
de chaleur sensible gagnée ou perdue dans l'unité 
de temps par l'unité de volume du fluide suivi dans 
son déplacement éventuel a dû passer par conduc- 
tibilité thermique à travers les surfaces qui la 
limitent. 


() C. R. Acad. Sc., 1941, 218, p. 988; 1942, 214, p. 301; 
1942, 245, séance du 3 novembre. 
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Re 


IV. 


On sait que par l'unité de surface perpendiculaire 
à un axe Ox passe une quantité de chaleur propor- 
tionnelle au coefficient de conductibilité À et au 


: : oT 
gradient de température = 


Découpons dans le fluide un élément de volume 
cubique ou prismatique mobile dont les arêtes, de 
longueur unité, sont parallèles à trois axes Ox, Oy, Oz. 

Si l’on considère ses deux faces opposées, séparées 
par l’unité de distance dans la direction Ox, il entre 
par la première, dans l’unité de temps, la quantité 


D 21 


- 01 c 
de chaleur À qui sort par la seconde, augmentée 
nus : 2EN 2 : 
de sa différentielle 1 par rapport à x. La chaleur 
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nette gagnée ou perdue entre les deux faces consi- 


o 


o? 
dérées est donc À A 
À ce flux calorifique gagné ou perdu dans la 
direction Ox s’ajoutent les flux gagnés ou perdus 
dans les directions Oy et Oz. 
La variation totale de chaleur par unité de volume 
dans l’unité de temps est donc 
COTE T ES OET, 
Ds de dy? T 92 } 


Cette variation de chaleur sensible se traduisant par 
une variation de température dT' dans le temps df, 


où p représente la 


elle est encore égale à co _ 
masse spécifique et c la chaleur spécifique par unité 
de masse. 

Égalant les deux formes de la chaleur sensible 
perdue par l’unité de volume du fluide dans l’unité 


de temps, nous obtenons 


Es CÉTESRO RE 22 
(jé dt = dx? dy + Dal . 


En général, le fluide est en DRE et l’on devra 


par conséquent calculer la dérivée, — L de sa tempé- 


rature en l’accompagnant dans son Re 
- La température T dépend alors non seulement du 
temps { comme en un point de coordonnées fixes 
données, mais aussi des coordonnées %, y, z à un 
instant donné de la masse considérée; sa variation 
totale en un point entraîné par le fluide est donc 

oT oT OT OT 


dE DA an He 


Divisons par df pour faire apparaître sa dérivée, en 


- 
notant que 7 ne sont autres que les compo- 


, 4, ds 
OTTAET 
santes U, DL, W Fe la vitesse du fluide. 

Ainsi apparaît l'équation classique de distribution 


des températures 


= mt 72 +P —+- (Y - 


0T' OT oT OT 
( dt dx dy ds ) 


(4) 


, RU AT 
ax? dy? fe 5) A 


Les termes dits de convection sont ceux où figurent 
explicitement les vitesses u, v, w; ils représentent la 
contribution que le mouvement du fluide apporte à 
la distribution des températures dans son sein. 

Cette équation reçoit en général une démonstration 
différente. Considérant un volume dx, dy, dz fixe 
dans l’espace, on écrit qu’à la chaleur d’échauffement 
près, le fluide y apporte par ses composantes de 
vitesse la chaleur qu’en emporte la conductibilité. 
Cette vue s’identifie évidemment à la précédente, mais 
sa présentation prépare moins bien à ce qui suit. 

Pratiquement, on ne considère que le régime per- 
manent, tel que la vitesse et la température du fluide 


soient stationnaires en chaque point. La différen- 


tielle "5 est alors nulle. 


. Le coefficient de convection. — La chaleur 
qui, en régime permanent, passe de la masse fluide 
à une paroi solide qu’elle baigne, est directement liée 
à la distribution des températures T au voisinage 
même de la paroi. 

En effet, cette quantité de chaleur doit transiter 
par conductibilité à travers la couche fluide qui 
reste immobile au contact immédiat de la paroi, 
quel que soit le type d'écoulement considéré. 

Si l’on nomme Oy l’axe normal à la paroi, orienté 
vers le fluide, la quantité de chaleur qui est ainsi 
transférée du fluide à la paroi (y — o) dans l’unité de 
temps et par unité de surface commune a pour 
expression 
2 OP: (5) . 

2 ‘OV 7 3=ù 

Il est commode de rapporter cette quantité de 
chaleur à une différence de température caractéri- 
sant le champ thermique étudié, par exemple à la 
différence T'.—0, qui existe entre la température T. 
prise loin dans le fluide et la température 0, prise 
à la surface du solide. On pose alors 


01 
OO ,, = = FRE ; 
=) = ae 0) 


Ceci suppose un fluide plus chaud que la paroi. 
Si le fluide est plus froid que la paroï, (TT. —6,) 
et () changent simultanément de signe, et les 
y=0 


dy 
équations précédentes représentent encore bien la 


chaleur ‘qui est transférée de la paroi au fluide par 
unités de surface, de temps et d’écart thermique, en 
l’admettant négative. 

Le facteur « ainsi introduit a été appelé « coefficient 
de convection ». Il caractérise l’aptitude du fluide à 
transférer la chaleur, indépendamment des diffé- 
rences de température, puisque le rapport 


ee (5) 
À sta. HE 6 dy y=0 


est sans dimensions par rapport aux températures. 
Si les températures T, et 0, varient d’un point à 
l’autre de la surface, il faudra évidemment envisager 
une définition locale de «. 
Pour un écoulement à une seule direction Ox, on aura 


7 (Tr). 
Le = ——— : 
; (Te — 60 x dy 0) 


La définition (2) du coefficient de convection peut 
être con$idérée comme absolument générale, tandis 
que l'équation de la chaleur (1) ne l’est pas. 


3. L’'équation généralisée de la chaleur. — La 


chaleur sensible cp … » qui dans l'unité de temps 
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apparaît ou disparaît dans l’unité de volume du 
fluide, peut bien y être venue ou l'avoir quittée par 
transfert à travers les surfaces-limites, mais elle 
peut aussi ÿ avoir été créée ou dissipée entre ces 
limites : si le fluide est le siège d’une réaction 
chimique ou d’une transformation physique qui 
y dégage ou qui y absorbe de la chaleur; s’il s’y 
produit une diffusion qui substitue du fluide plus 
ou moins chaud à du fluide plus ou moins froid; 
s'il est soumis à un rayonnement qui le chauffe ou 
le refroidit directement, etc. 

Appelons d’une façon générale + Qq la quantité 
de chaleur qu’apportent ou qu’empruntent à l’unité 
de volume et par unité de temps toutes ces nouvelles 
causes énumérées, et toutes autres que l’on pourrait 
imaginer. 


L’équation complétée de la chaleur s’écrit évidem- 


ment 
JT Te CAM o0T ; D LMO IRON: er œ 
dt dx TR en 3 D ne T ) 


q étant généralement une fonction de *, y, z, U, v, w 
2 X À 

Le double signe + rappelle qu’en plus du cas, 
essentiel en pratique, où une réaction exother- 
mique tend à élever la température du fluide par 
un g positif souvent très sensible, il faut envisager 
aussi les cas des réactions endothermiques, des diffu- 
sions froides, des pertes par rayonnement, ou même 
des évolutions physiques vers des états moins con- 
densés (vaporisations, fusions), qui tendent à l’abais- 
ser par un g négatif. 


4. Formes possibles de la chaleur sensible 
créée ou détruite. — Explicitons q dans quelques 


cas, avec les approximations permises en pareille 
matière : 


19 Le rayonnement réciproque entre une couche 
fluide à température T et des parois opaques à 
température 0,,< 1 comporte une double infinité 
de réflexions et d’absorptions. Si le fluide est iso- 
therme et n’offre qu’un constituant doué d’absorp- 
tion, avec une « épaisseur réduite » / (produit de 
sa concentration partielle + par la pression totale p 
et l'épaisseur réelle L) (?), il y apparaît, par unité de 
surface parallèle aux parois et par unité de temps, 
une quantité de chaleur rayonnée : 


Re 
De CE Pt un 


(?) La notion d'épaisseur réduite postale que les degrés 
d'activation des gaz sont autonomes, et qu’en particulier la 
présence de gaz « inertes » n’influe pas sur le degré d’acti- 
vation des gaz émetteurs, ce qui n’est pas toujours stric- 
tement vérifié. à 

(3) Voir M. VÉRON, Tendances actuelles des techniques de ia 
chaleur, p. 34 et suivantes. 
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Es, et E7 sont les pouvoirs émissifs du corps noir 


aux températures 0, et T': 
UE AR gr-cal 


T+ à 
) m?>x h x vKt 


avec s— 4,91 
l 


as, est le pouvoir absorbant (ou émissif relatif) 
des parois; si on les assimile à des corps noirs : 


CAES 


0 


a,,est le pouvoir absorbant de la veine sous l’épais- 
seur réduite /, à la température T : 


> “ Con D JE dÀ 
bandes ne nee TELE Ê À 
pi EF d'A 


UT rss 
Le signe DA rappelle que l’émission des fluides 
bandes | 

procède par bandes de longueur d’onde, générale- 
ment très étroites, dans chacune de ces bandes, la 
valeur moyenne de a,, dépend très peu de T, en 
sorte qu’on peut sortir de l'intégration le facteur où 
elle intervient. 

On a donc assez sensiblement : 


rpeS De Ii un)! | ———— 


bandes 


à 
hs 


= = He) 


la fonction f (!) étant une fonction croissante amortie 
de !, et o(T) une fonction généralement décroissante 
de T, qu’en première approximation on écrira T-$, 
avec £ > o. 

Si plusieurs constituants rayonnants sont méêlés 
sans réagir, leurs pouvoirs absorbants a, s'ajoutent. 

D'autre part, le volume gazeux mis en cause par 
unité de:surface d’une couche d’épaisseur réelle L 


est égal à L = - 


La chaleur qu’apporte (ou qu'emporte) le rayon- 
nement, par unités de surface et de temps, peut 
donc s’écrire : 


MAC ; 
RE D 


fluides 


OT) 


où F(l) — aie est une fonction décroissante de 1. 


En fait, la convection empêche la veine d’être iso- 
therme, mais l'expression (4) s'applique avec une 
approximation souvent suffisante à chacun des élé- 
ments de filets pris avec sa température T propre. 


2° La diffusion d’un autre fluide, plus chaud ou 
moins chaud, se traiterait comme il suit : 


Soit 7 le « taux » du nouveau fluide par rapport 
au fluide principal, et soit D le coefficient de diffu- 
sion. À travers l'unité de surface perpendiculaire à Ox 
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me RUE OT! 
et dans l’unité de temps passent par définition D sp 


unités de volume du nouveau fluide. 
L’aprort de ce fluide dans l’unité de volume par 
unité de temps sera donc 


227 227 d2% 
“ ( Der À 
Si, en se substituant au fluide principal, l’unité 


de volume du nouveau fluide apporte + q calories, 


il Vient 
92% DE (3 
Top Lo 02? Se 


d27 
<q gb (DE + 


Le transfert de chaleur par turbulence, qui s’appa- 
_ rente au cas précédent, sera étudié plus loin dans 
un chapitre spécial. 


39 Pour le cas, très important en pratique puisque 
c’est celui des flammes, où la chaleur naît dans le 
fluide du fait d’une réaction chimique, nous raison- 
nerons sur un exemple simple mais facile à géné- 
raliser. 

Considérons un écoulement dans la direction Ox, 
parallèlement à ‘une paroi, et admettons que des 
molécules d’un type unique y agissent sur d’autres 
molécules d’un type unique pour donner des molé- 
_ cules d’un composé défini unique, par exemple 
Hi+=(Os+4N:) = H:0+92N. 

Soit n le nombre des molécules primitives d’un 
type (H? par exemple), qui sont encore présentes 
dans l’unité de volume du fluide au temps f et à 
l’abscisse x, sur le filet de vitesse u; soit Æ q, leur 
chaleur de réaction moléculaire. 
dn ; 
se molé- 
cules (« vitesse de réaction »); elles « actualisent » 


Dans l’unité de temps se combinent — 


= de lories 
—- TE Ca 1esS. 

On a évidemment 

du 
TI=+g: EF 

Il faut connaître l’évolution de n en fonction de f. 

Nous admettrons que la vitesse de réaction varie 
par exemple suivant une loi 


Lio 


a étant une fonction supposée connue de la tempé- 
rature T, qui dépend elle-même du temps t. 
On en tire 
L 
— a(t) dé 
qu = — a(t) dt, i j 


= — puis An —= No € 

mn représentant le nombre des molécules initiales 
qui intéressent l'unité de volume à l’instant quel- 
conque f. Il est égal au nombre N, des molécules 


LU 


initialement comprises dans cette unité de volume, 
raréfié dans le rapport inverse de la dilatation subie 
par le volume initial. 

Pour les liquides, la dilatation est négligeable ; on a 


10 = Ne 


Pour les gaz, on peut souvent admettre qu'en 
raison du ballast d’azote apporté par l’air, la réaction 
est équimoléculaire et qu’ainsi le volume actuel suit 
la loi des gaz parfaits; on a 

SE 
no= No a . 

{ est le temps pour un observateur qui accom- 
pagne les molécules dans leur mouvement; suivant 
la direction x, sa différentielle s’écrit donc 


QC = ee 


Le nombre des molécules actives encore présentes 
dans l’unité de volume à l’instant considéré apparaît 
ainsi égal à 


1, = fa on. 
IL — ae T € 


et leur chaleur de réaction RES la forme 


De Led 


DE GhES SE GO SE qaN À 


Le groupement g.N, est la chaleur potentielle 
totale des molécules qui sont initialement contenues 
dans l’unité de volume; c’est une donnée essentielle 
du fluide considéré; nous la désignerons par Q; > 0, 
puisque le signe convenable figure explicitement 
dans la formuk précédente. 

Finalement, la chaleur dégagée par unités de volume 
et de temps prend pour expression 


(6) 


: T à 
où le facteur m0 valable pour les gaz, est à Supprimer 


pour les liquides... 

Cette expression de q à une valeur très générale et 
elle s’applique à toute loi de réaction propre aux 
molécules considérées, car la température T et le 
nombre actuel n des molécules étant liés par l’évo- 


, lution de la réaction, la loi que nous avons admise 


pour la vitesse de réaction équivaut à une loi d’aspeet 
plus général 


et elle englobe donc toute loi particulière des ciné- 
tistes, telle que 
dn 


= 47 — b1P n1. 


L 
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5. Position du problème concernant le coeffi- 
cient de convection. — Par les expressions (4), 
(5), (6) de q, nous formons des équations de la 
chaleur (3) qui généralisent l'équation classique (1). 

Elles ne sont peut-être pas nouvelles en elles-mêmes, 
et il est à souhaiter que les termes introduits aient 
déjà été considérés effectivement toutes les fois qu’ils 
existaient, encore que nous n’ayons vu figurer le 
« terme complémentaire » q que dans l’équation de la 
conduction accompagnée d'effet Joule, conduction 
qu’on pourrait aussi qualifier de vive. 

Mais, ce qui est essentiel, les conséquences con- 

cernant le coefficient de convection défini par l’équa- 
tion (2) n’en ont pas été tirées. A notre connaissance 
tout au moins, leS physiciens n’ont calculé que des 
coefficients de convection répondant par la. for- 
mule (2) à des distributions de température régies 
par l'équation restreinte (1). 
_ Leurs expressions, largement diffusées, et les valeurs 
numériques qu’en déduisent les ingénieurs, supposent 
ainsi réellement, sans qu’on y prenne garde, que la 
chaleur transmise vient du fluide préalablement 
chauffé et est apportée uniquement par son mouve- 
ment, pour être transportée finalement au solide par 
la conductibilité du fluide. 


Or, d’après la relation (2) qui les définit, les coefñ- 


cients de convection sont conditionnés par les gra- 
dients de température à la paroi, qui eux-mêmes sont 
régis par le «terme complémentaire » g, concuremment 
avec les facteurs classiques À, €, p, U, v, w, &, z. 

Toutes choses égales par ailleurs, il existe donc 
autant de coefficients de convection qu’il y a de dis- 
tributions possibles pour la chaleur q créée par unités 
de volume et de temps, c’est-à-dire qu’il y a de 
sources de cette chaleur, et même de combinaisons 
de ces sources. 

En particulier, le coefficient de convection d’un 
fluide dépend de son rayonnement, ce qui proscrit 
déjà en principe l’existence d’une expression ‘uni- 
verselle du coefficient de convection, valable aussi 
bien pour un fluide pratiquement non émissif 
(comme l’air) que pour un fluide nettement émissif 
(comme une fumée). 

Inversement, le rayonnement d’une veine fluide 
dépend en principe de sa distribution de tempé- 
ratures, et par conséquent de ses échanges par 
convection. 

La conception habituelle qui ajoute les flux de 
chaleur échangés par convection et par rayonnement 
calculés indépendamment l’un de l’autre est donc 
fausse, mais elle est trop commode pour qu’on 
l’abandonne. 

Les couches fluides immédiatement voisines des 
parois offrent d’ailleurs. avec elle des écarts ther- 
miques, restreints (nuls à la limite pour y — o), de 
sorte qu’en fait leur rayonnement réciproque avec 
ces parois est limité, n’influençant donc qu’assez peu 


la distribution locale des températures au contact, . Au du 


qui régit le coefficient de convection. 
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Nous nous proposons de déterminer d’abord les 
coefficients de convection déduits de l'équation 
généralisée de la chaleur (3) dans le cas de l’écoule- 
ment laminaire, avec le terme complémentaire q 
fixé par l’expression (6), c’est-à-dire pour un fluide 
en réaction, le rayonnement et la diffusion étant 
encore négligés; ces formes moins massives de 
la chaleur créée feront l’objet d’une étude ulté- 
rieure. 

Afin de rendre notre exposé plus intelligible en 
facilitant les comparaisons, nous avons jugé utile 
de résumer auparavant le travail classique fondé 
sur l'emploi de l’équation restreinte (1). 

Nous avons évité d’y introduire d’entrée la notion 
de « couche-limite ». On verra que si l’on s’en tient à 
ce qui se passe dans son épaisseur, la forme du 
coefficient de convection n’est pas affectée; seules 
ses constantes sont légèrement modifiées, en convec- 
tion vive comme en convection morte; mais il était 
désirable de n’en pas préjuger. 


IL. La convection morte d’un îiluide passif 
s’écoulant en régime laminaire permanent, 
le long d’une plaque. 


Nous rappelons que nous entendons par fluide 
« passif » ou «mort », un fluide plus ou moins chaud, 
dont la température peut varier, mais dans lequel 
la chaleur ne peut entrer ou sortir d’un élément de 
volume mobile que par le mécanisme de la conduc- 
tibilité. 

Dans un tel fluide, l'équation de la chaleur est 
l’équation restreinte (1), dont la résolution suppose 
connue la distribution des vitesses. 


1. Solution de Blasius pour la distribution 
des vitesses. — Soient Ox un axe pris dans le sens 
de l’écoulement le long d’une surface cylindrique ou 
assimilable, et Oy un axe perpendiculaire à sa 
surface. 

P représentant la pression, & la viscosité absolue 
et W la fonction potentielle des vitesses, les équations 
classiques de l’hydrodynamique des fluides visqueux 
incompressibles s’écrivent d’une façon générale 


ou du 9P 5 Au du 
» SLT Le Sn (a) 


I 
7 + RE = 
dx dy 0 0x 2 Vox? dy? 
de dep 1 OP u /@p 6 
U—— HP = —: He — + — J) (db) 
dx dy 0 0Y o \ dx? dy? 
ou de è do , CL e 
= —0, AVEC Un — ee 
dr dy 1 dy ë dx ù 
On a évidemment 
PEU, 


dre < dy? ? 
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x qu 
ce qui permettra de négliger les termes où entrent > 


et les différentielles de v (4). 
D'après l’équation (b), où celles-ci figurent à tous 


; à )P 
les termes sauf un, le gradient de pression normal =. 


devient négligeable; la pression P n’est plus Re 
que de x; elle se trouve donc imposée par un phé- 
nomène extérieur au fluide. 

L'hypothèse fondamentale formulée par Blasius 
est qu’à grande distance de la paroi, les forces de 
frottement dues à la viscosité p n’agissent pas; la 
vitesse u ne dépendant plus de y devient là u.(æ), 
et ses différentielles par rapport à y disparaissent. 

L’équation (b) donne alors 


d’où, en substituant dans les équations générales à 
distance finie, 


Pr PANNE RE ER («&) 
ox dy dx e dy? ? ù 
NON Dhs EN CE (c) 
07 ? dx 
AVEC OM OUR OI POU TETE 0 


Puisque nous avons supposé (en c) le fluide incom- 
pressible, nous devons logiquement admettre qu’il 
ne se dilate pas (liquide ou gaz presque isotherme), 
ce qui, pour un écoulement plan, entraîne évidemment 


ou 2 


He (tt) Const: donc 7 
OT 


0): 


Il reste alors pour (a) l'expression fondamentale 


du du CE ne 


= y 


“oz ‘dy Ip ? 


(a) 


en posant 


H 


IV 


(viscosité cinématique). 
P 


Considérons, avec Blasius, les deux facteurs : et & 
tels que 


UNI Ce (7) 


he NE Ve et 


On constatera que & n’est fonction que de £, ce qui 
nous permettra js décrire le champ de vitesses avec 


la seule variable 
Vz 
Il faut expliciter £(E). 


(‘) On pourrait même négliger v, comme nous le ferons 
dans l’étude thermique. 


On doit avoir d’après (b) 


D'une part 


© 
fx 
|| 

[Di 

4 

A 

œ 

œ 

A 
Il 
Te 
£ EM) 

8 
" 
ES 
& 

Il 
"nr 
R 

s 


soit 
ue 
U= — C: 
Dr et 
(8) 
I “Us ñ 
D / (EE —C) 
cp} (A 
Par conséquent 
5 
du Loto ye ef 6 
Vas Tr) A ETES 
dr » 5 2 
COTE or J ce 
HN — 4: où 6%; —= = , 
dy PAU à aV vx 
0? [22 (LÆ* Fo 12 QT) s 2 
et CL A + f! x où EU 0 


2 y? D: 4 


Et l'équation hydrodynamique (a') s'écrit 


u? YU VU re 
Le GLCEER + UE (he yet, = detre, 
4 
, d’une part, 
o AA 
u=, Us YU 
ee AVE LR LCTE, 
4 CCE ï — €: Æ 
2 2 
us / ny us # " 
! = USE CREUSER IEEE 
5 GAE RAGE 


et, d'autre part, 


2 » 

(4, VU (US STD EE u? 
er = CRE, = > CHR et TA 
‘ A #4 * T 5 ns = 4x 


Il reste donc simplement 


ter =— tr, (a) 


avec les conditions aux limites suivantes : 


L — (Age 
ue \Y =0, donc A0) 
rec, donc ER, 
, (u—r—=0, donc Ciel, donc os 
il faut ; 
l Abe donc Cr 


La solution &(£) satisfaisant à ces conditions s’ex- 
prime par un développement en série 


Cat 
Ca OR EE PES 
nr Êe Er) A d 
U=0 (9) 
où 
CREME ECS EN TANT ES 37, Ci 97806 


CRU 
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Sa dérivée première 
N=N 


Car , 
— 220, Nm ARR EE) 10 
u= > 1) (3n +1) IR Fe de 


PE 


permet de calculer les vitesses u et v par (8). 

Il faut prendre &« —1,326 pour réaliser %: — 2 
pour £ — , 

Notons encore en passant que la force de frot- 
tement f par unité de surface, sur la face balayée 
de la plaque, est 


du Ur pr 


avec Den et FR 
= = F 
; \ 


HIS 0 (11) 


Ainsi, sous la condition que le fluide ne se dilate 


pas et qu’en conséquence sa vitesse u, à l'infini 
reste uniforme, la distribution des composantes de 
vitesse u et v, régies par (8), ne dépend des coor- 
données x et y que par la fonction *(°), qu’exprime (9), 


_ de la variable auxiliaire : que définit (7), c’est-à- 


dire qu’en définitive les vitesses s'expriment bien par 
l’unique variable 


F= 


= Z. (12) 


Al 


\ 


que nous retiendrons pour la suite immédiate. 
En d’autres termes, les courbes 


Us g y 
— — = const., ou Z= “= — const. 


I 
2 ADN Tr V L 


é(æ, F)= 


qui sont des paraboles ayant leur sommet sur le 
bord d’attaque, figurent des lignes d’égales vitesses u 
et v, ce qui rend sensible la signification physique 
de € (fig. 1). 

Si l’on suppose que le frottement est localisé au 
voisinage de la plaque, dans une « couche-limite » au 
delà de laquelle se conserverait la vitesse de l’écou- 
lement potentiel, on trouve pour son épaisseur, en y 
admettant la distribution de vitesse u(y) suivant 


Blasius : 
= 6T RE 2 6/2 C ie 
(US 


D’après (7), elle se trouve coïncider avec la 
courbe £ = 3. D'après (10) qui donne € et (8) qui 
donne # cette épaisseur correspond en réalité à 


= ,Gs=1-e 
| 


avec € très faible mais non nul. 


(5) Au lieu de la constante 6, on a parfois adopté les con- 
stantes 5,45 ou 4,92, qui répondent la première à un pro- 
fil u (y) simplifié, l’autre à un choix arbitraire de &. 


Partant de zéro pour x = 0, 9, atteint prati-. 
quement, aux échelles usuelles de x, quelques milli- 
mètres pour les liquides et quelques centimétres 
pour les gaz. 


À 
l FLUIDE A L'INFINI 


Distribution des vitesses dans un écoulement laminaire 
le long d’une plaque, suivant Blasius. 


Lignes d’égale vitesse (4 — const.), couche-limite (£ < 3). 


La remarquable propriété qu'ont £ ou Z de régler u 
et v cesse d’être exacte dans tous les cas où uw, 
varie (sauf un qui n'offre pas d'application réelle); 
mais l'expérience a permis d'admettre qu’elle repré- 
sente encore une approximation pratiquement satis- 
faisante dans les cas effectivement rencontrés. 


2. Solution de Pohlhausen pour la distribu- 
tion des températures et le coefficient de convec- 
tion laminaire morte. — Pohlhausen observe que, 
comme celle des vitesses u et v, la distribution des 
températures T' ne doit dépendre des coordonnées x 
et y que par la variable unique 


He 


Z= — 
\æ 


(41) 


En effet, en régime permanent le long d’une plaque, 
on peut poser 
ZX MIE UE PT PT PT 
NN ae LÉ: = à ù 5 g 


et l'équation de la chaleur (1) devient 


($) La suppression de cette hypothèse compliquerait les 
calculs sans changer leur nature et sans modifier apprécia- 
blement les résultats. 
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Posons que T' ne dépend que de la variable Z. 
On doit avoir 


Î 


oT 70 7 PT 


ER E PAR DE 


et, x s’éliminant, l'équation de la chaleur prend la 
forme 


I (0 HR s 
SCPUZ SEAT 0: (43) 


Or, les facteurs c, 0, À sont des constantes ou assimi- 
lables à des constantes, et la vitesse u atteinte à une 
abscisse x donnée ne dépend que de la variable Z. 

Introduisons les deux grandeurs sans dimensions 


“e À [42 
S — — (nombre de Stanton) ET (A) 
o) Us 


n ne dépend comme u que de Z, si u, ne dépend 
pas de x. 
L’équation de la chaleur devient alors 


TEL Ad di 7 


dZ? —= 0, 
et elle fait bien apparaître T comme une fonction de 
la seule variable Z. 

Elle équivaut à 


(4) 

AE 

dZz US PRE 

HET NCAA 
az 


Immédiate, l'intégration donne 


COUPÉS (Z)Z 47 

73 — — 12) 47 
d7 no : 

7 = AeE 
dZz 


donc évidemment 


CS 
a) 4 


On en déduit aussitôt le gradient des températures 
dans le fluide à la paroi 


Gin G), E)- 4 
FÈe dy 10 dz 2-0 Vx ; 


qui, transporté dans la formule de définition (2) 
du coefficient de convection, donne l'expression 


PEL) ue NPA 
ET 18% Ed 04 dy Fi) E Væ T, — 0 x 


On détermine la constante d'intégration À, et 


même directement le rapport — en exprimant 
come 0: 

qu'intégrée depuis la plaque (y — 0) jusqu’au loin 

dans le fluide (y — «), au droit de l’abscisse x, la 

différentielle de température dT doit restituer l'écart 


thermique total T, — 0, relatif à cette abscisse. 


PHYSIQUE N° 5. 


Il vient 


LES 
IT J EN 
in PA) dz 


del 


JE 


ht) "ie 
VS 
— A [ e 0 
7==0 


2 
n(2)2 dZ 
dZ =7T,— 0, 


d’où 


ne | 
rs (| PAPA 
Do dZz 


124 F, 


Le coefficient de convection s’explicite donc comme 
il suit 


où les facteurs Z, dZ et n(Z) s’exprimerit tous en 
fonction de la seule variable auxiliaire : de Blasius. 
D'après les relations (7) et (12) comparées, on a 


donc CAT (15) 
LIL EE AS 
es 


= 
AT) "(© et 
US ; 


D'après les relations (14) et (8) combinées, on a 


tu ll ‘ 
N — == — £S. (16) 
44 2 + 


En remplaçant dans l'expression précédente du 


coefficient de convection, il vient, après simplifi- 
cations, 


soit encore 


(17) 


avec 


AE (18) 


où l’on remplacera &: d’après (10), ce qui donnera 


Su 


5 5 Casa 
CAE — ( IVe du-+2 QE 
ie CAE : > ) Gas is 


0 


soit 


$ » È Canxi Can-+S 
Lt 1)" É A € 
J OR. ) (3 +1)! 3 +3 EEE 


No 5. 


La formule (17), complétée par (18) compte tenu 
de (19), est l’expression exacte du coefficient de 
convection d’après Pohlhausen. 


Dans une Note antérieure basée sur la consi- 


dération de la « couche-limite », nous avions ét écon- 
duits à l'expression 


don ue 
Ar — a 
ANT 22 ANT A 


4 
a —T [°n(0)6 40 
VF = | e sJ d8 (5) 


0 


(17) 


avec 
(18°) 


où la variable auxiliaire O0 représentait la cote y 
rapportée à l'épaisseur 5, de la « couche-limite » 
à la même abscisse x 


Z 0 Y 1 Le 
(] = —— = —, <— 
Ve at 
64/ 


alors qu'ici 


Cette variable auxiliaire Ü était donc liée à la 
variable = considérée ici, par une loi de proportion- 
nalité 


ainsi qu’on peut le voir en partant de la significa- 
tion physique de €, l'épaisseur de la couche-limite 
correspondant à £ — 3 et à Ur. 
En fonction de :, il vient 
u ne + dE £ d£ 


A = Ù = -= 0 d0 — 
ï] sn 2 d D d0 : 


SAUT di CE dE; 


ce qui conduit à 


2 # "= 


Pour, 
(2 


D'où, en négligeant l'influence du fluide au delà 
de la couche-limite, entre 0 — 1 et 0 — 0, 


61, DUT 


ce qui entraîne l’équivalence des deux expressions (17) 
et (17') du coefficient de convection 
Il faut calculer l’intégrale purement numérique J;, 


définie par la formule Se fe t:c d: étant fournie 


par la formule (19), l’on donnera à n, successi- 
vement, les valeurs o, 1, 2, etc. 


(?) C. R. Acad. Sc., 1942, 21%, p. 301. 

Pour éviter toute confusion avec la variable de Blasius, 
à laquelle il est très généralement affecté, nous remplaçons ici 
par le symbole 6 le symbole © utilisé dans notre Note en cause. 
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On obtient 


Fe . 
11,426)? 66 11(1,826)% 61 
re 
ONG FAT Gi 


Mettons en évidence le groupement variable 


1 1 
1,326 \7 38 \5 
e 2 F 5 nl 
g ( 39 ) LL ARR (5) de 


On obtient 


Le 7 2 2 
3 S 2 aP—= 2 g + = CLS SP 
dE e 2.4! 7! ds. 
1,326 " 6 
@ (] 


Lorsque la variable est petite, les termes en #6, 
0%, ..., qui entrent dans l’exponentielle sont négli- 
geables devant le terme en v$. Lorsque ® est grand, 
l’exponentielle elle-même est négligeable. 

En consentant à une erreur secondaire qui sera 


corrigée plus loin, on se bornera donc à écrire 


Ji = | ——— S emtido,. 
7: (x) Ê jh ne 
1 


3 3 2 
(5) 1,319 et ji e de 0,852, 


0 


avec 


d’où, enfin, 
1 | l 
JMD 107802 ST L20H OS 2/1 © 9,941 8". 
En définitive, le coefficient de convection admet 
pour expressions 
à : ny 
VA Us RTE te [12 
Lp= _— — 0,408 LP X == 0 


2,2418 


En faisant disparaître la viscosité cinématique » 
au profit du nombre de Stanton $ 


J /. 


L À . 7 57 Acpus 
sn soc 


D » d’où 


on peut encore lui donner les formes suivantes 


! 
“S6 |COU 4) À u 
= A/R ossi / RE. (21) 
DTA D æ 


Le calcul complet de J,, compte tenu des termes 
supérieurs en +, a fourni à Pohlhausen un coeffi- 
cient 0,332 au lieu du coefficient 0,446 : c’est là 
l'effet des termes négligés. 

Bien que ces termes fassent intervenir $ un peu 
différemment, nous admettrons que les formes (20) 
et (21) restent valables. 

Les expressions plus exactes et classiques (5) du 


(#) Voir M. G. RisauDp, Mémorial des Sciences physiques, 
XL, p. 26 à 28. 
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coefficient de convection morte sont ainsi les sui- 
vantes : 


au de AUX ce 
Un 01932 0 00PX HE 0,332 S° Kepie + (22) 
væ æ | 


Notons qu’elles répondent à 


1 a. 1 

2  o,446 te 

JC TNIPONSEX 44 1,508", 
0,332 


donc à 


En particulier, pour les gaz biatomiques, qui 
dominent dans l’air et dans les fumées, on a 


LA 
cree 


Fa 9; 907; 


== 


Sr, LS 0,3328 180,301, 


et le coefflcient de convection se précise ainsi 


3. Garactère local des solutions de Blasius 
et de Pohlhausen dans les gaz soumis à des 
échanges calorifiques. — La solution de Blasius 


i 
KL (2) (1) 
à NN SZ 
._ REGION OÙ SE RSR 
DÉFINIT L'ÉCOULEMEN NN 4 
TANGENT RC N RQ 
KR 
TA D A ail te 
MP A ET NS 0 _ 
Fig. 2. — Définition d’un écoulement de Blasius 


tangent à l’écoulement réel pour une tranche d’abscisse Ax. 


pour les vitesses u et celle de Pohlhausen pour les 
température T supposent essentiellement que la 
vitesse « à l’infini » u, est constante et uniforme, 
c'est-à-dire qu’à quelque distance de la paroi elle 
est indépendante de la longueur + parcourue par 
le fluide le long de la plaque considérée. 

En revanche, la solution de Pohlhausen ne suppose 
pas explicitement, en plus, que la température T., 
à quelque distance de la paroi est constante et uni- 
forme. En effet, l'équation (2) qui définit le coefficient 
de convection «. peut s'entendre localement, pour 


No 5, 


une tranche Ax découpée dans le fluide autour 
de l’abscisse x envisagée, la valeur de T. et celle 
de 0, se rapportant à cette abscisse ; ceci, à la seule 
2e DÉTeRT Tee RE 

condition que 5 soit bien négligeable, comme nous 


2 


T 


avons dû le supposer. 


Pour que la vitesse à l’infini u, soit uniforme sans 
que la température à l'infini le soit, il est nécessaire 
que le fluide ne subisse aucune dilatation ou contrac- 
tion, ce qui n’est exactement vérifié que dans les 
liquides et les gaz isothermes, sièges de réactions 
équimoléculaires. 

C’est là une restriction fort importante, puisqu’en 
principe elle exclut des applications les cas les plus 
fréquents où le fluide, gazeux, cède ou perd beaucoup 
de chaleur. s 

Elle semble avoir échappé aux usagers, qui 
appliquent les formules (21) à (23), aux fluides dila- 
tables comme aux autres. 

Nous pensons que la solution de Pohlhausen reste 
cependant parfaitement suffisante pour définir un 
coefficient de convection à l’abscisse x, si l’on intro- 
duit dans les formules (21) à (23) la valeur numé- 
rique de T. et celle de u, qui conviennent à la 
situation locale de l'écoulement. 

Cette façon de procéder revient (fig. 2) à remplacer 
l'écoulement réel, qui, dans son ensemble, diffère 
de l’écoulement de Blasius à cause des dilatations 
ou contractions, par une suite d’écoulements de 
Blasius dont chacun offre localement un profil u(y) 
de vitesse (courbe 2) qui aboutit, sur la paroi, tangen- 
tiellement au profil u(y) réel (courbe 1). 

Réalisant le même gradient de vitesse (5) à la 

JA 
paroi, ils respectent chacun le frottement local 


du 
f=t ie 


totalisant ainsi le même frottement d'ensemble. 
D’après l’équation fondamentale de la chaleur, 

ils respectent aussi chacun le gradient de tempé- 

OT 


rature ( ) à la paroi, donc le transfert de chaleur 
Y—=N 


22 
: OT 
DEMI 
ù : ( dy jo 


unitaire local 
totalisant ainsi le même échange d'ensemble. 

Pratiquement, nous appellerons écoulement de 
Blasius tangent relatif à l’abscisse x celui qui, sur 
une tranche Ax aussi mince qu’on le voudra prise 
autour de +, donne le même frottement local que l’écou- 
lement vrai; ce sont ses caractéristiques, notamment 
sa vitesse à l'infini u,, connue par le frottement 
grâce à la formule (11), que nous aurons à introduire 
dans le calcul du coefficient de convection. 

On profitera de ce fractionnement pour ajuster 
les valeurs de S$, c, p, À ou p dans la couche-limite, 
qu'on a dû supposer constantes bien qu’elles varient 


A 


Î 
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avec la température, sauf » pour les He El et S 
pour les gaz. 

En résumé, la solution de Pohlhausen pour le 
coefficient de convection x. n’a de valeur absolue qu’à 
la limite, pour des fluides non dilatables; dans les 
autres cas, il n'a qu’une valeur locale, valable pour 
une tranche d’abscisse Ax d'autant moins étendue que 
l'isothermie est moins bien vérifiée, encore que l'usage 
se soit établi de l’étendre dans tous les cas à l’ensemble 
de la surface balayée. 


III. — La convection vive d’un fluide actif 
s'écoulant en régime laminaire permanent 
le long d’une plaque. 


1. Possibilité d’une solution pour la distri- 
bution des températures. — [L’équation de la 
chaleur devient (3) au lieu de (1), la distribution 
des températures T en fonction des coordonnées x 
et y est donc autre. Pour la même température à 
l'infini, le profil T(y) des températures à l’abscisse x 
(fig. 3)se déforme de (1) en (2), sa pente est modifiée, 

s À oT ; 
en particulier la valeur (Tr). qu'elle prend sur 
la paroi. Le coefficient de convection «., qui est 
toujours défini par la relation (2), n’obéit donc plus 
à la même loi. 

On se rend tout de suite compte de l'erreur qui 
consiste à appliquer les formules (21) à (23) à la 
transmission de la chaleur par une flamme ou par 
une lessive en réaction. 

Pour un fluide, soit chaud (T. > 0), soit froid 


| (T < 05), le profil des températures (fig. 3) devenant 


soit (2), soit (2) au lieu de (1), le gradient de tempé- 


: 14 À 
rature sur la paroi, » et par conséquent le 
à =} 


( 
72 
L 1 7 ÔT 
coefficient de convection, &,—= =—— | — > se 
1 00 \9Y /3— 


trouvent être augmentés ou diminués, suivant que 
la réaction échauffe ou refroidit le fluide au voisi- 
nage de la plaque, pour une même De 14 
à l'infini. 

Nous nous proposons essentiellement de calculer 
un coefficient de convection qui puisse toujours se 
comparer directement à celui de Pohlhausen. 

Il semble dès l’abord qu'une difficulté insurmon- 
table s’y oppose : 

Dans la solution de Pohlhausen, bâtie sur l’équa- 
tion restreinte de la chaleur (1), la température T 
ne dépend que de la seule variable Z — Je mais 

V2 4 
quand, dans l'équation généralisée (3), nous intro- 
duüisons le terme complémentaire q, qui dépend sépa- 
rément de x et de y, cette circonstance simplificatrice 
disparaît. Pour déterminer le problème, il faudra se 
fixer des conditions aux limites différentes, plus pré- 
cises, de sorte que la comparaison envisagée avec le 
coefficient de Pohlhausen perdra toute signification, 


à supposer que l’on ait pu intégrer l'équation de la- 
chaleur. 

Écrivons l'équation générale de la chaleur (3) 
dans le cas de l'écoulement plan, en exprimant que, 


Fig. 3. — Distribution des températures à l’abscisse x 
dans le fluide, pour une température à l’infini T< donnée : 


Profil (1) en convection morte, 
profil (2) en convection vive (== q étant du signe de To — 65). 


Les deux coefficients de convection sont entre eux comme 
les pentes, sur l’axe normal y, des tangentes aux deux profils 
sur la paroi. 


pour les termes renouvelés de l’équation restreinte (1), 
la température T ne dépend que de la variable 
unique Z. Il vient, par analogie avec (13), 


PP NE OC 0): (24) 

Si nous admettons que l'écoulement de Blasius 
se conserve, la vitesse u n’est fonction que de Z, 
mais le terme du deuxième membre, æq(u, x), n’est 
pas réductible à cette seule variable. 

Cependant, nous avons vu que, si l’on y regarde 
de près, le coefficient de convection de Pohlhausen 
n’a pour les gaz qu'une signification locale, propre 
à une abscisse x déterminée. 

Si l’on consent à raisonner encore localement, 
pour les abscisses d’une petite tranche Ax autour 
de x, il est possible d'exprimer le deuxième membre 
de (24) en fonction de la seule variable Z, après 
quelques transformations tout à fait acceptables. 

La connaissance de l’écoulement de Blasius, jointe 
à celle de la nouvelle fonction æg (u, x), permet alors de 


aT 0 
procéder au nouveau calcul de +» donc de ( Fe) Ù 
7—0 


et d'en déduire, par une Fnéthéde calquée sur celle 
de Pohlhausen, un coefficient de convection «; 
qui est localement exact et directement comparable 
à celui de Pohlhausen pour la même abscisse x, 
comme correspondant au même type d'écoulement 


et défini en fonction des mêmes éléments. 


2. Explicitation du terme complémentaire 
dans l'équation de la chaleur. — Nous consi- 
dérerons le cas particulier très important où la 


2 


chaleur qQ apportée à l'unité de volume provient 
d'une réaction telle que la combustion dans une 
flamme, et s'exprime comme on l’a vu par (6), Q; 
représentant ici la chaleur potentielle contenue (ou 
encore contenue) dans l’unité de volume du fluide 
quand il aborde la plaque sur l’arête x — 0. Le terme 
complémentaire admet donc pour expression brute 


1 


= / PAVA) de 


; ; À 
rq(u, 2) = Qraxr D e © 


Son aspect suggère de substituer à la variable x 


a dx 
[A me 


. . . AL . È 
les variables sans dimensions y? Puis : 


, ii a “) ax “ 
a autre 
[7 21 == { [l = 
1 (D 0 U (ke | cost [7 


CT AN O p Cp e ‘0 
a dx 


en écrivant 


On posera 


DS 


NM= — (20h) 


et 


À 4 dx 
ni dz y: / / « : 
Po — CEA = jé (26) 
0 É ù 


Il vient alors 


721 


TON AIONPT ce I M. 


Dans chaque tranche Ax autour de l’abscisse x, 
la vitesse u varie rapidement à partir de zéro et la 
température absolue T varie moins rapidement à 
partir de 6,, quand le filet considéré s’élève à partir 
de y — 0. 

Les groupements m, et m, varient évidemment 
peu, ce. qui nous autorise à les remplacer, dans la 
tranche considérée, par leurs valeurs moyennes 
obtenues en faisant varier y. Nous grouperons en une 
seule constante locale l’ensemble des termes «moyen- 
nisés » 


IR — y M. (2) 


Si la température ne varie pas considérablement 
d’un point à un autre, le facteur a(T') varie assez peu 
lui-même, et la quantité m, définie par (25) reste 
très voisine de 1. 

Nulle pour x = 0 et pour 4 = æ, la quantité m, 
définie par (26) passe pour y = 1 par un maximum 
T 
2,718 

ce Maximum. 

C’est seulement pour les valeurs très faibles ou 
très fortes de 


le DIE (8) 


0 u U 


k ns ee 
égal à = - Elle reste longtemps voisine de 
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c’est-à-dire pour les réactions presque achevées(n ->0), 
les vitesses de réaction très lentes ou très rapides 


(4 -> o où x), les vitesses d'écoulement très rapides 


ou très lentes (u- > © ou o), les abscisses balayées très 
courtes ou très longues (x-> 0 ou +), que m, devient 


‘ 11 
très petit par rapport à son maximum ? - 


Le long d’une plaque donnée, x croissant de zéro 
à l’infini, n décroissant de n, à o, la constante locale m 


définie par (25), (26) et (27) part donc de zéro pour 


y revenir, après avoir plafonné sous son maximum, 
sur une Jongueur Ÿ Az d’autant plus étendue que la 
vitesse de réaction est plus faible et la vitesse 
d'écoulement plus forte. 

D'importantes erreurs dans l’appréciation de la 
loi de réaction ne peuvent modifier beaucoup les 


composantes m, et m, du facteur m, donc ce facteur 
lui-même, ce qui est satisfaisant. 

Pour les cas normaux de flammes dans les fours, 
on ne se trompera guère en posant en première 
approximation 


_— —— 1 I 


NE = MIN: —= 3 ou ve 


et il sera toujours loisible d'obtenir une valeur 
locale plus précise de m par des relevés de tempé- 
rature et des analyses de gaz pris dans les flammes. 

En définitive, dans chaque tranche Ax, nous assimi- 
lerons le terme complémentaire à sa valeur moyenne, 
qui s’écrira 


SEULE) = Om Er (2) 


7 


Bien entendu, nous déformons un peu le profil T(y), 


rt 


SES (2 
d’où une certaine erreur sur (7) » donc sur la 
à y—=0 


valeur du coefficient de convection qui sera calculée 
plus loin. 


3. Solution généralisée pour la distribution 
des températures et le coefficient de convection 
laminaire vive. — [L’équation de la chaleur (24) 
prend localement la forme 

PEU ENRRS AT LH D: 
VAT de 122000 Qy Co PA 


DA 


Pour les liquides, le facteur .; 


est à remplacer 


par r et o — ps, densité initiale. 
Pour les gaz, le facteur à 
de :, et alors o est variable. 

En appelant p, la densité normale et p, la densité 
initiale des gaz de la flamme, qu’on assimilera à des 
gaz parfaits, il vient 


est en général différent 


" 
Lu 


at” 


D 
Î 


LES 


11e 


/ 


N° 5. 
donc 
1 7 CAN 1 nt Ti m 
LR = —=— DRE 0 EU = ee Op 4 
6 T Pn Go e + ce JD u QE cher 
Introduisons la variable auxiliaire = définie par (15), 
telle que 
" LE: % 
AE) VE ë 
soit 
aol d7' de r7. 427 
de En de oi de 


ainsi que les deux grandeurs sans dimensions définies 
par (14) 

ve % u 

O = — et = —": 
CE ls 


L’équation de la chaleur devient, pour tous les 
fluides, 
Î mQr 
+ LR Ra) U, 
SAR) 


(30) 


où ?, désigne la densité initiale pour les gaz comme 
pour les liquides. 

Ainsi qu’on l’a vu, n ne dépend comme u que de :, 
si u, ne dépend pas de x. On a encore 


1,, 
= — Es. 
DS 


L'’équation de la chaleur (30) fait donc bien appa- 
raître T comme une fonction de la seule variable €. 


Son intégration est possible, et la solution générale 
CT) 
dE est 


d’où évidemment 


CYA 
Fe = ne 


On en déduit aussitôt le gradient des températures 
dans le fluide à la paroi 


es 2 (9 (%) =1/#0 
()..= D , dE E—0 en VIT 


qui, transporté dans la formule de définition (2) du 
coefficient de convection, donne l’expression 


à 7 RD ES Dre D) 
LATE dl’ ma 0o Le je 53 » Ve V4 TA, FL CP 
On détermine la constante d'intégration D, et 


; D : 
même directement le rapport 7,5, exprimant 


qu'intégrée depuis la plaque (y = 0, : — o) jusqu’au 


loin dans le fluide (y —æ, : =), au droit de 
l’abscisse x, la différentielle de température dT doit 
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restituer l'écart thermique total T, — 5, relatif à 


* cette abscisse. Il vient 


a = Ë=» 3 = ee te 6 4ë 
ue ïE _ E dE = D Je sf dé 
EN /—0 


> mQr [—* ET ; 
Poe “A dé e me 
o Cpo t—0 
— LE 0o, 
d’où : 
LA 1 Sole at 
ÿ RS et OT fr) 
1% C 0(L'» eo 0) E—0 À 
{ m2 2 
E 2 f'utar 
| CRE 5) ae | 
D + 0 " i 
RTS É 


Le coefficient de convection s’explicite 
comme il suit 
À Phi Dane mOr a 
An = —— — es —_— 1 
de Tee 5 Re ce 
en posant comme antérieurement 
= NT 
7e [ e So al (48) 
CA) 
et, d'autre part, 
PANNE ALNCEe à Sid 
JL 1 de ee tee d£|+ (32) 
20 


pe S | 


La formule (31), complétée par (18) et (32) compte 
tenu de (19), est l’expression localement exacte du 
coefficient de convection généralisé. 

Dans la Note antérieure des Comples rendus, 


basée sur la considération de la couche-limite, 
nous avions été conduits à l’expression 


à SET ALOr, | : 
= —— — D 1 
“ av PTE He PME 
avec, comme on a Vu, 
1 
62e 3S" 
et, d’autre part, 
; (] 
DOUTE ) ü ds 0 1 f (66 a 
0 dôe 5 de n(B)eS* dé, (32) 
0 


où 


ue  1(0)9 d6 — LE dE, 


[< Edé 


donc 


d 
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_ce qui conduit à 


one A 1 d : 
I 1 — À He dù ee < A b:S 45 
D del NO ie 
= 0 


(0 


En négligeant l'influence du fluide au delà de la 
couche-limite, entre 0 — 1 et 0 — +, il vient 


% 
È 


ce qui entraîne l’équivalencé des deux expres- 
sions (31) et (31') du coefficient de convection |. 

I faut calculer l'intégrale purement numérique J; 
définie par.la formule (32), 1e C6 d£ étant fournie 


par la formule (19), où l’on donnera à n, successi- 
vement, les valeurs o, 1, 2, etc. 

Comme pour l'intégrale J,, la solution exacte est 
très laborieuse. 

Compte tenu des autres approximations faites, il 
est indiqué de s’en tenir au premier terme du dévelop- 
pement de &, d'où 


% 326, 
fear = 2e, 
© 


= Apr CL TE 
Ja = | dte *%$ 1 1,326E e °® : 
E0 0 


Mettons encore en évidence le groupement variable 
3S 


1 
1,326 : RS DT 
es E, avec dé — (5) do. 


On obtient après simplifications 


ce qui donne 


0 = 
T 


=» ? 
2=3s f dpe-® f à 
Et (2 


La valeur de J, est ici proportionnelle à &. 
Pour le dormbre S — 1,33 des gaz biatomiques, on 


e?' do. 


S 
a très sensiblement ——— et des calculs très 


1, 
1,326 
fastidieux, sur lesquels nous ne nous étendrons pas, 


donnent 
JP I, 3% 


Pour l’ensemble des fluides, nous pouvons donc. 


admettre comme valeur approchée 


Vs 1,69 


Nr 


100 


En définitive, le coefficient de convection prend 
pour expression 


(00) 


RTBE mor à 
Se l ma 1 10) | [ ay [ax]» 


On remarquera que $ n'apparaît plus au second 
terme du crochet. 
En faisant disparaître la viscosité cinématique » 


No à. 


au profit du nombre de Stanton $, on peut encore 
lui donner la forme suivante 


k 
Li 0, Hs, Re 2129 30 


Pour les gaz biatomiques en particulier, le coeffi- 
cient de convection d’après (33) se précise ainsi 


2 mOf 
Ay= 0, 3OLA VE nr STE Te AN) 


Les calculs qui nous ont conduits à ces formes du 
coefficient de convection vive n’impliquent aucune 
hypothèse restrictive concernant l’évolution possible 
de la réaction, antérieurement à l’arrivée du fluide 
sur la surface (x — o) : seule, la constante Q; doit 
être éventuellement réduite de la chaleur déjà libérée 
par la réaction. 


mr 
coo(T'>—00) 


4. Caractère local de notre solution géné- 
ralisée. — Comme la solution de Pohlhausen con- 
cernant |[«./,, et pour les mêmes raisons, notre solu- 
tion concernant «. suppose que le fluide n’est pas 
sujet à 
n’est exactement vérifié que dans les liquides et dans 
les gaz isothermes. 

Quand le fluide est le siège d’une réaction, telle 
qu’une combustion, elle accentue généralement les 
variations de température à distance dans le fluide, 
encore qu'il existe des cas où elle les atténue au 
contraire. 

Comme la solution de Pohlhausen appliquée aux 
fluides dilatables, notre solution reste cependant 
suffisante pour définir le coefficient de convection 
local à l’abscisse x, si l’on introduit dans les expres- 
sions (33) ou (34) les valeurs numériques des fac- 
teurs $, À, €, p, u. qui Œœnviennent à la situation 
locale de l’écoulement, et notamment la vitesse à 
l'infini u, d’un écoulement de Blasius fangent à 
l’écoulement vrai pour cette abscisse. 

Le critère pratique d’un tel écoulement sera encore 
qu'il donne la même perte de pression que l’écou- 
lement réel, sur une tranche d’abscisse Ax, comme 
il a été dit. | 

Si ce critère est lui-même inaccessible, on pourra en 
première approximation définir l'écoulement tangent 
en lui donnant la vitesse et la température à l'infini 


de l’écoulement réel. 


Par la suite, dans les cas où l’on pourra craindre 
qu'à distance de la paroi, l’écoulement réel s’écarte 
trop de l’écoulement tangent, on pourra demander 
à des expériences systématiques les coefficients de 
majoration/ou de minoration qu'il convient d’ap- 


_pliquer aux facteurs u, et T. de l’écoulement réel 


pour obtenir des résultats corrects. 
Mais nous devons observer que ceux qui jusqu'ici 


E CARRE 


| = fs + [or a | 


des dilatations ou à des contractions, ce qui. 


PS 


(3) 


34) 
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ont appliqué la théorie de Pohlhausen, et a fortiori 
toutes autres plus simples, n’ont pas connu de tels 
scrupules. 


En résumé, avec ou sans réaction, la solution de 
Blasius pour l'écoulement n’est valable que dans 
un domaine restreint, et la vraie signification des 
formules (33) à (35) est de permettre, en chaque 
abscisse, une comparaison locale légitime de la con- 
veclion vive ou totale avec la convection morte de 
Pohlhausen, ce qui constitue un résultat essentiel. 


9. Expressions approchées de l'écart ther- 
mique total. — En principe, pour le cas pur d’une 
convection avec réaction, l'écart thermique total 
T.—%% de l'écoulement réel résulte de l'équation (3) 
de la chaleur écrite au delà de la couche-limite, q 
restant exprimée par (6); on a 


rs 
te 
0 


; avec 4 = Qyan pe e 


L'intégration pour T', exige l’introduction d’hy- 
pothèses plus grossières que celles ayant présidé au 


. Calcul de «, par ( œ)- = 


Pour les liquides, =" ne joue pas,ona u.—u,—const. 


< . 
et p. — Po = Const. 
Dans la mesure où a. est invariable, il vient 


ce qui donne 


Dia (es 
f1£s — Bo = (76 — 065) = Qt + 


He 
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: 1 
Pour les gaz, 7 joue, u., p. et a. 


sont très 
variables. 

Si l’on admet que la dilatation des filets se fait 
uniquement suivant Ox, et que a.(T.) est linéaire, 


hypothèses courantes des théories classiques de la défla- 
gration, on a 


1h To £ 
A Fe = 07 Lot lies 
donc 
de CN 
— = — — const 
LUE üo 
et il vient 


ce qui donne 


en 
PERTE) RS fee ï). (200 
C Pa 

Sinon, il faut procéder par tranches. 

En fait, la réaction donne souvent naissance à des 
constituants qui, comme CO? et H?20 dans les 
flammes, sont nettement émissifs et entretiennent 
avec les parois un rayonnement réciproque qui atté- 
nue fortement les variations de l’écart thermique à 
quelque distance, surtout si l’écoulement est lent ou 
la réaction « paresseuse »; la diffusion de veines 
froides et les dissociations jouent dans le même sens. 

On n’est donc presque jamais en droit d’expliciter q 
à l'infini, simplement d’après (6). Nous reviendrons 
sur les calculs plus exacts. 

Pour l'instant, notons que dans le cas des flammes 
de fours, les observations montrent que l'écart T.—0 
reste souvent sensiblement uniforme-sur des surfaces 
étendues, sous les influences combinées de la com- 
bustion et des dissociations, du rayonnement et des 
diffusions. 

(A suivre.) 


GROUPEMENT ET DÉGROUPEMENT AU SEIN D'UN FAISCEAU CATHODIQUE 
INJECTÉ DANS UN ESPACE EXEMPT DE CHAMPS EXTÉRIEURS, 
APRÈS AVOIR ÉTÉ MODULÉ DANS SA VITESSE. I. 


Par R. WARNECKE, J. BERNIER et P. GUÉNARD. 


Sommaire. — {Lorsque l’on provoque par un moyen quelconque, dans un faisceau électronique, une 
perturbation périodique du temps de transit, il se produit un groupement spatial des électrons; la densité 
de ce groupement dépend, en particulier, de la forme de l’action perturbatrice. Dès que ce groupement 
apparaît, il se manifeste, dans le faisceau, des forces de répulsion d’origine électrostatique qui tendent 
a modifier le rassemblement en voie de formation. 

Dans la première partie de cette étude, les auteurs considèrent des tubes électroniques dans lesquels 
un faisceau rectiligne, après avoir été modulé dans sa vitesse, est injecté dans un espace où toute force 
électrique d’origine extérieure est éliminée grâce à une électrode étui formant écran, et ils déterminent, 
sous l’hypothèse du rassemblement cinématique, la forme de la modulation de vitesse à adopter pour 
obtenir un groupement conduisant à une cession aussi complète que possible de l’énergie cinétique 
des électrons à un champ haute fréquence sinusoïdalement variable agissant sur les électrons paral- 
lèlement à leur mouvement à la sortie de l’électrode étui. 

Dans la seconde partie du travail, les auteurs étudient, pour le même type de tube, la forme que 
revêt théoriquement le groupement quand on tient compte des répulsions interélectroniques. Pour simplifier 
les développements mathématiques, l’analyse est restreinte au cas d’un faisceau cylindrique parfai- 
tement focalisé de l’extérieur, ayant un diamètre faible devant celui de l’électrode étui et subissant une 
modulation de vitesse sinusoïdale de faible profondeur. 


Introduction. 


0. À. — Aperçu des problèmes. — Comme il 
ressort de plusieurs exposés [1] [2] de l’un des 
auteurs, les principes fondamentaux des principaux 
tubes à modulation de vitesse récemment développés 


pour la production, l’amplification et la détection 
des oscillations électromagnétiques d’ultra-haute 
fréquence peuvent être établis à partir du schéma 
de la figure 1 correspondant aux tubes dits à « con- 
version par glissement ». Un faisceau rectiligne 
d'électrons est injecté suivant l’axe d’un cylindre, 
les électrons étant accélérés avant leur entrée dans 
ce cylindre et collectés après l’avoir quitté. Le canon 
électronique et le collecteur d'électrons n’inter- 
viennent pas dans les effets de haute fréquence 
produits, et il est inutile ici d’en envisager la forme 
Constructive. 


Le fonctionnement de ces tubes repose sur les 
actions suivantes : 

Aux bornes d’une section MM' d’un tuyau conduc- 
teur (organe modulateur ou « rassembleur »), on 
maintient une différence de potentiel périodique 
qui, alternativement, accélère et ralentit par rapport 
à leur vitesse moyenne, les électrons qui traversent 
cette section. Cette action fait apparaître une pulsa- 
tion de densité dans le faisceau, c’est-à-dire une 
modulation du courant de conduction qu’il trans- 
porte. Ensuite, par induction, la composante alter- 
native de ce courant, au niveau de la section CC’, 
cède de l'énergie à une impédance Z connectée aux 
bornes de cette section (organe collecteur d’énergie). 

Dans la théorie la plus simple de ce système, 
on établit les lois du groupement des électrons dans 
l’espace séparant MM’ et CC’ en négligeant les effets 
des électrons les uns sur les autres, et, pour sim- 
plifier encore, on suppose qu’il n’existe dans cet 
espace aucun champ électrique alternatif ou constant 
d’origine externe, si bien que le groupement se fait 
par « glissement ». Une telle théorie, appelée 
«théorie du rassemblement cinématique », esquissée 
par Heil [3], développée particulièrement par 
Webster [4] et les auteurs [5] [6], est complétée 
dans la première partie de cet exposé. 

Cette théorie qui tire son intérêt de la simplicité 
de ses bases et de sa claire représentation physique, 
est évidemment insuffisante pour rendre compte des 
faits expérimentaux puisqu'elle néglige tous les effets 


sù 
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de la charge d'espace et, notamment, ceux qui 
résultent de la naissance des « paquets d'électrons » 
dans le faisceau; pour conduire à une image plus 
correcte du phénomène, elle doit donc être au 
moins complétée ou corrigée : c'est ce qui est fait 
dans la seconde partie de cet exposé. Après avoir 


Calculé par approximations successives les forces 


qui s’exercent sur les électrons et les modifications 
consécutives de leur mouvement, on détermine, pour 
un certain nombre de cas pratiquement importants, 
l’amplitude et la phase des premiers harmoniques de 
l’onde du courant de conduction transporté par le 
faisceau. 

Dans l’une et l’autre des deux parties dé l'exposé 
on supposera implicitement : 


a. Que les champs électriques agissant sur les 
électrons à l'extérieur de l’espace de groupement 
sont développés entre des électrodes infiniment 
voisines l’une de l’autre. De cette façon, on néglige, 
a priori, dans les espaces correspondants, les angles 
de transit des électrons, dont les -auteurs ont par 
ailleurs [5] [6] déterminé, sous les hypothèses du 
rassemblement cinématique, l'influence particu- 
lièrement importante dans le collecteur. 


b. Que les tensions alternatives créant les champs 


qui agissent sur le faisceau sont appliquées entre 


des électrodes idéales dont les surfaces, perpendicu- 
laires à l'axe du faisceau, sont équipotentielles. 
L'introduction de cette restriction est importante : 
non seulement elle fixe la direction des forces élec- 
triques qui s’exercent sur les électrons, mais elle 
élimine encore la nécessité d'introduire, dans les 
calculs, des coefficients d'efficacité des électrodes 
tenant compte de ce que leur action sur un électron 
dépend de sa distance à l’axe du tube. 


0. B. — Expression générale de l'intensité du 
courant de conduction transporté par un faisceau 
dans lequel le temps de transit est modulé. — 
Considérons un faisceau rectiligne et « homogène » 
d'électrons se déplaçant dans le vide. Soient 4, 
l’époque de passage à l’origine MM et { l’époque 
en un point P. La quantité d'électricité tdi, qui 
est passée en MM’ entre les époques & et & + dé, 
passe en P entre les époques t et { +- dt; d’après le 
principe de conservation de l'électricité, on a : 


A) dé —= it dé 
soit 
Û I J 
RE | 
2) dé d7' ( ) 
dt dt 


où T —{—1, est le temps de transit des électrons. 
Si celui-ci subit une perturbation périodique, c’est- 


à-dire si T est fonction périodique de f, (de pulsation 
fondamentale w), : sera aussi fonction périodique 
0 
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de {, (donc de f) et pourra être considéré comme une 
somme d'ondes sinusoïdales 


avec 


EURE 1 DT d(wt). (2) 


a ; 
— R(Æneirwt), 
27 : Lo 


Je I 
Lo D 
En tenant compte de (1) et en faisant le change- 
ment de variables d'intégration £— {, + T(to), 
avec ol, —9, o=+r,ona 


kn PTE e—julg+tig) de. (3) 
: F 


Les k, s’obtenant par intégration de (3), le calcul 
de à revient à celui de 7 — r(+) en tenant compte 
des facteurs déterminants qui, dans le cas général, 
ont une origine à la fois extérieure etintérieure au 
faisceau. 

Ce résultat, tout à fait général, qui met en relief 
l'effet du temps de transit, d’une part est indépendant 
de la forme des trajectoires et de la forme du 
courant i, et; d'autre part, ne suppose rien ni sur la 
cause de la variation du temps de transit, ni sur la 
forme de celle-ci (pourvu cependant qu’elle soit une 
fonction périodique de t). 


PREMIÈRE PARTIE. 


Rassemblement cinématique des électrons 
dans un faisceau 
après modulation dans sa vitesse. 
Groupement optimum. 


1. A. — Cas d’un faisceau rectiligne injecté 
dans un espace exempt de champ de îorces 
extérieures, après avoir subi une modulation 
sinusoïdale de sa vitesse. Rendement de conver- 
sion maximum correspondant. — Quand on ne 
tient pas compte de la charge d'espace, le calcul 
du temps de transit est immédiat : en effet, soit 
v, = Vy(i + sine) la vitesse des électrons à la 
sortie de MM' (fig. 2); le temps que met un électron 
passant en MM" à l’époque t{, pour aller en P 


est TT; ()— * d’où, en introduisant les nota- 


= U Z 
tions na =013, 1= — 
(L 


î à 
PÉTER. 

Généralement, pour simplifier le calcul (et aussi 
parce que ce cas est pratiquement avantageux), on 


admet que 92à est négligeable devant 97, alors 


uæ À — Ô) sino, (4) 
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ce qui donne, d’après (2) et (3), 


— = — ein(g+k-üising) ds = entr J,(nôh), (3) 


re RS NE 

rie 2 rt RE) cos n(ot— À). (6) 
pl 

[Dans cet article, J,(x) représente la fonction de 

Bessel de première espèce, d'ordre n, et d’argu- 


ment x.] 


A 


a NE 4 
L 
‘ 
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Fig. 2. 


On appelle rendement de conversion le rapport 
de l’énergie cédée au champ ultra-haute fréquence 
du collecteur à l'énergie continue nécessaire pour 
accélérer le faisceau. 


Lorsque le champ du collecteur est sinusoïda] 


et que le faisceau est modulé en densité par le méca- 
nisme ‘ci-dessus, il résulte de la relation (6), ainsi 
qu’il est bien connu [4] [1], que la valeur maximum 
du rendement de conversion ne peut être égale qu’au 
maximum de J;(9), soit 0,582 (obtenue pour 


dÀ = 1,84). 
1. B. — Groupement amélioré grâce à une 
modulation de vitesse non sinusoïdale. -— Pour 


dépasser la valeur précédente du rendement de 
conversion, il est nécessaire d'abandonner la forme 
* sinusoïdale de la modulation de vitesse ('). La forme 
convenable peut être trouvée directement à partir 
de l’expression générale du rendement de conversion. 

Si dans la section CC’, d’abscisse ! à partir du 


() Ce problème de l’amélioration du groupement dans un 
faisceau commandé par modulation de vitesse, qui vient 
de faire l’objet d’une communication importante de BorGnis 
et LEDINEGG [9], a été traité antérieurement à ce travail 
par les auteurs du présent exposé; on trouvera, en effet, 
le résumé de la première partie de notre article dans le 
préambule théorique d’un brevet français déposé le 19 sep- 
tembre 1942 sous le n° 471.736. 


No ÿ. 


modulateur (fig. 2), le faisceau électronique, dans 
lequel la modulation de vitesse a provoqué une modu- 
lation de densité, traverse un champ électrique étroit 
dont l'intégrale linéaire dans le sens du mouvement 
a pour valeur instantanée —U, cos wt, le rendement 
de conversion pour le réglage optimum de phase 
a pour expression (°) 


TJ céleste (7) 
"0 


où. U, désigne la tension continue correspondant 
à la vitesse moyenne v, des électrons, et 7, leur 
angle de parcours entre le modulateur et le collecteur. 

En général, la valeur maximum qu’il est possible 
de donner à U,est U,, car, pour une valeur plus 
élevée il y a réflexion d’électrons dans le champ 


(3) L'expression (7) peut être établie indépendamment de 
la forme du courant transporté par le faisceau de la façon 
suivante (voir [r}) : Soit à, l’intensité instantanée du courant 
transporté par le faisceau au niveau du plan médian du 
collecteur. Si les électrodes du collecteur sont très voisines 
l’une de l’autre, le courant induit qui traverse l’impédance 
shunt de ce dernier est sensiblement égal à à,. D’autre part, 
la tension instantanée que ce courant développe entre les 


A de 
électrodes est u, — U,cos (wf + B)et, d’après (1),ona &=to FT 


Il s’ensuit que la puissance cédée par le faisceau au champ 
électrique est 


27 
W = iu,= = f LU, L° cos(wt + B) d (ue). 
0 


On obtient l’expression (7) en divisant W par i,U,, en rendant 
nul l’angle de phase B et en remplaçant wf, par ® et wl 
par @ + t (+). 

Au sujet de l’expression précédente de la puissance utile, 
il convent de noter : 


a. Que le fait d’aimettre une form: sinusoïdale pour x, est 
légitimé si l’on suppose que cette tension est développé: aux 
b:rnes d’un sy.tème ayant les propriétés d’un circuit anti- 
résonnant. 

b. Que la composante du courant induit due au dépla- 
cement des éle:trons est li1 seule qui présente un intérêt réel 
puisque la composante cpacitive du courant total est 
éliminée dans le circuit extérieur et ne produit aucune puis- 
sance moyenne (voir [2], chapitre II). 

c. Que pour une durée non négligeable du transit des élec- 
trons dans le collecteur, une expression de la puissance utile 
pourrait être trouvée et utilisée comme celle qui est donnée 
plus haut. La valeur instantanée du courant induit à intro- 
duire pour son calcul (valeur du courant induit résultant des 
courants induits élémentaires dus à tous les électrons présents 
à l’instant considéré { entre les électrodes) s’obtiendrait à 
partir de l’expression 


l 
‘dx 
LE À UA : pi(i) 9; (£,) ‘dt dé,, 


œu ((4)) 


dans laquelle 6,, v,, 4 dési nent respectivement la densité, la 
vitesse et l’époque de passage des électrons au niveau de la 
première électrode du collecteur et T la durée du transit 
(fonction du temps) entre les électrodes de celui-ci, distantes 
de h; maïs le calcul de la puissance sur ces bases conduirait 
alors à des résultats d’une forne si compliquée que leur dis- 


‘cussion présenterait difficilement un intérêt pour le problème 


considéré ici et c’est pourquoi le cas du temps de transit 
négliseable ‘dans le collecteur est seul considéré. 


Î 


No 


du collecteur et perte correspondante d’énergie. 
D'autre part la valeur de l'intégrale est maximum 
lorsqu'on a, quel que soit 0: 


cos[ S + s/(9)] =, soit Et) A7. 
Ainsi, théoriquement, dans l’hypothèse du rassem- 
blement cinématique, le rendement de conversion 
d’un tube à temps de transit peut êtré amené à une 
valeur voisine de 100 pour 100 (quel que soit le 
mode de groupement) en faisant en sorte que +, (œ) 
soit une fonction périodique de + coïncidant 
avec 2kK7— dans une période. Ce résultat est 
général, il est indépendant de la forme de la trajec- 
toire des électrons et du mode de groupement. 
Dans le cas particulier du groupement par glis- 


2 GW ns 
sement, on. a 7; — et la condition de rendement 
1 


de conversion maximum est évidemment, en 
* w L 
posant À:= — 
Po 
Po 
o+— =92kr. (8) 
P1 
1. C. — Forme de « tension » de commande 


rendant maximum le rendement de conversion 
dans le cas du groupement par glissement. — 
Si l’on continue à supposer que la modulation de 
vitesse dans le faisceau est provoquée par une diffé- 
rence de potentiel alternative appliquée entre deux 
grilles idéales très voisines l’une de l’autre (couche 
éléctrique double infiniment mince) portées au 
potentiel U, constant, la vitesse v, de l’électron à 


.la sortie du système de commande correspond à 


2 
un potentiel U, + U(o) et l’on a (2) = I + _. 
0 0 
La condition de rendement maximum s'écrit alors, 
dans le cas du groupement par glissement : 
IE (=) —: 
Un : 
Soit 


do + 1 COS O + D1 Sin? +...+ anCosn 9 + br sinnp+.….., 
\ 


le développement en série de Fourier de (=) ; 
1? 

4r2K(K—1) 

valeur moyenne de U(o) étant nulle, a, doit être 

égal à l’unité; ceci impose à À, une valeur particu- 

lière À4 


le calcul donne a, — + D'autre part la 


U,= 27 VA(K— 1). (9) 
Les coefficients a, et b, (cf. Annexe I) sont donnés 
par 

2} 
472 k(k — 1) CA 


br PM Foi(a(k— nn) — Citaknr)], 


nn 


[sé( 2(k—1)nr) — si(2knx)], 


dn = 
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où si et Ci sont les fonctions sinus-intégral et 
cosinus-intégral. Ces coefficients tendent vers zéro 
lorsque n augmente indéfiniment et, en tenant 
compte de (9), on peut les représenter par les for- 
mules asymptotiques 


DR DE ADN 
An TT IE ET) 
PES ut e 
—92 x 2! DEN 
bar —— ? PRE 
my, n3} 7 
Application numérique. — Pour un tube fonc- 


tionnant avec U, — 3 000 V sur une longueur 
d'onde de À — 25 cm, la longueur de l’espace de 
glissement [e = \2 VE (k — D) serait de ! — 6,6 cm 


pour k—=3 (À,—= 15,4 radians), de: { = 9,35 cm 
pour kK—4 (1,— 21,8 radians), de [= 12,05 cm 
pour k—5 (l, — 26,2 radians); dans ce dernier 
cas les coefficients du développernent seraient 


& À 0,015 
As & 0,004 


bn —0,142/n, 


An © 0. 


Pour obtenir le rendement maximum de conver- 
3 UE 
Sion n = jy l'intégrale du champ dans le collecteur 
à 


étant toujours de la forme U, cos wf, il faudrait 
donc que la tension de commande dans le modu- 
lateur soit représentée par 


ù =(— 0,140 sing + 0,016 cos?) 
: + (—0,071 sin2® + 0,004 cos2%) 
4 0,142 . 
— 0.047 Sin 39 —..:.— Smart 
n 
1. D. — Forme de la vitesse à la sortie du 


rassembleur rendant le rendement de conver- 
sion maximum dans le cas du groupement 
par glissement. — La vitesse , doit satisfaire 
à l'égalité (8) et doit avoir comme valeur moyenne v, 
(fig. 3). 
Posons 
ï | 
— = 1+ di COS® 
Po 


+ b!, sins 


+ a, cosnp + D, sine +... 


: PE , (9 sp 
La condition que.la valeur moyenne de " soit égale 


à 1 entraîne pour À, la valeur À; déterminée par 
Lau ! 

= f kr dy, 

2 Rate ou DAT CESR 


SE: 
= M [logo kx — log2(k —1)r], 


c’est-à-dire 
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d’où 1. E. — Forme de l'angle de transit 102) 
RIT RENTE d rendant maximum le rendement de conversion 
Fe ART AN dans le cas d’un mode de groupement quel- 

og 


a | 


Les coefficients a, et b, sont donnés par 


I Fe COS NN Pb cosn 
"nc < t °0$ 
ii = f À Ans Rae IE #1 do, 
= 2kKr— 9 7 Ù 
0 T (A —1)R 1 
ST <1 S(Â—1)T à à 
pl sr Sinno Tac SEE de Nes 
ii A la Free CEE a 
RCA TC MARS RE É 
d’où 
ZT d’. 
= = Ci(oknz)— Ci(o(k—i)nz), 
VE 
r ; . 
= — sù (2 ANT) — st (2(A—i)uz). 
AYA É 
: U è 
Fig. 3. — Forme de rendant le rendement de conversion 


0 
maximum pour le groupement par glissement. 


De la même façon que précédemment on peut 
obtenir des valeurs approchées de a! et b’, à partir 
des développements Ai MEN ÈCE dé st ete CL, 
soient 


ra, I 2k— 1 
HORS A 1 : 2 I ee 
EU Mr ED) On) (en datlh 


Application numérique. Reprenant l’exemple 
précédent pour k — 5, on trouve comme valeurs 
des coefficients (3) 


0,071 0,000 


) = — =: 
n° D n° 


(5) Les a, et b’, ont des valeurs moitié de celles trouvées 


(®) 


dans le développement de Ü parce que k a été choisi 
0 


pour que le rapport des amplitudes des tensions soit petit; 
par ailleurs, les conditions de rendement maximum sont 
évidemment les mêmes, À, = À,, car 


_—- & 27 WA(K 3). 


conque. — On a déjà montré que la condition de 
rendement maximum, quel que soit le mode de 


groupement, est que la fonction périodique z/(0) 
de période 27, coïncide avec 2 kr — + entre o 
CUS, 


On en déduit le développement en série de Fourier 
de la fonction </,(9) optimum 


AR “ VAE 
ie) = a+ ai éoss + bT sine +...+ ascosns +b;sinne+.…. 


avec 
I PT 
= = He rc) do Orne, 
DR Je 
ja, I L 1: 
D co 
ed 
_— 
TMS 6 2 
= 1 (2k7 — ©) sin no do = —-; 
= far # 
01) 
d’où 


rC)—=(2k 1)" 


: ‘sin2v 
Er de SIN £ + 


Sin /à © : 
» 1 à 


1. F. — Rendement de conversion correspon- 
dant à quelques formes particulières de grou- 
pement. — Il est évident qu’au temps de transit 
optimum exprimé en angle par la relation (10), corres- 
pond, lorsque U, — U,, un rendement maximum de 
conversion égal à l’unité 


Er 
I ; co 
COS | 2+—2sinc+ Sing 
27, 4 : à 


Dre DS 
+ 3 SIN32+...+— sine. [de= 
È ' 2) : 3 


je 
1. (HOUis) 


Examinons ce qu’il advient quand le temps de 
transit n’a pas cette forme idéale, maïs une forme 
qui correspond seulement à l’utilisation des premiers 
termes du développement. 


a. Si le temps de transit varie sinusoïdalement 
[comme c’est le cas d’après (4) dans les tubes à modu- 
lation de vitesse à groupement par glissement, où 
la tension haute fréquence de commande est faible 


et sinusoïdale] le réglage défini par 
T9) =(2k—1)r +2sins 
donne, pour ,U, =, le rendement de conversion 


CS 


js cos (o-+ 2 sin?) do = J(2) = 0,577. 


I 


TU — 
DT 


Cette valeur est à peine inférieure à celle du 
maximum de r (0,582) trouvée au paragraphe 1. A. 


re 


NP=S 


et’ correspondant à un réglage optimum tel que (4) 
(2) =(2k —1)x +1,84 sino. 
b. Si le temps de transit varie comme asino+b sin2o 


le réglage défini par les premiers termes du dévelop- 
pement (10), soit 
(S)=(2k — ir +2 Sing + sine (12) 


correspond pour U, = U, au rendement de conver- 
sion. 


1 : - : 
COS (o + 2 Sinv + sin29) do 


edge 


D 
= — f [cos(o + 2 sin9)cos(sin2) 
— sin(o + 2 sino) sin(sin2o)] do 


Es s 3 = 
COS(S + 2 Sin fc ) + 2 dau (1)00s 4 ve | de 


Le sh 
CAT E 
LE t, 
Jill 
é Sr ) à 
Pire ss SIN(S + 2 Sin 9 [: Jap) sin rene | 
7 44 
| 


D 


on obtient finalement 


2 = | Jo(t)J1(2) —Y 240) Jig 02) Jig1(2)] 


1 


+Y op —1)Jop 101) Jp 2(2) 


XL 
et, en négligeant les termes inférieurs à 1074, 
= DhG)h(2)— 2086) +2%02(0) = 0,73. 


On peut constater que le maximum de l’inté 
grale | 


I 5 : 
cos(o + a sino + b sin2o) de 


ne diffère pas sensiblement de 0,737. Ainsi, à partir 
de deux termes, le réglage optimum est prati- 
quement défini par les premiers termes du dévelop- 
pement de (10). 


c. Si le temps de transit varie comme 


a sino + b sin2o +csin3v, 


le réglage optimum est tel que 


He) & RD +0 (sine + IE + PE) (43) 


(#) Lorsque le temps de transit dans l’espace de glissement 
ne correspond qu’à un nombre limité des termes du dévelop- 
pement (10), il convient, pour obtenir la valeur maximum 
possible du rendement de conversion (puisque des termes 
manquent alors), de donner aux amplitudes des termes restants 
des valeurs différentes de celles qui seraient adoptées si le 
déve'oppement complet était utilisé. 


(11) 
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et cela correspond pour U,= U, à un rendement 
de conversion maximum 


= 
I 
Use 


27% 


5 nee TA 
COS (e +2? Em" in3z) do, 
0 


dont la valeur approximative est 


M3 —= A3) à.) ou io) JA (3) 


x Jo(D[Z(2)—Ji(2)] 
+ (OZ 6)—-D()—J(2)] | 
+ BOB) (02)+ 70120) 0)1! 


+ Ja (:) JO)[B (0) — Ji (2)] & 0,767. 


Les résultats indiqués sous a, b, € montrent que 
le fait de produire, dans un tube électronique du 
type considéré, une variation du temps de transit 
de la forme (12) au lieu d’une variation sinusoïdale 
de la forme (11) apporterait, pour les conditions 
de réglage optimum, un accroissement du rendement 
de conversion de l’ordre de 16 pour 100 tandis que 
la forme (13) n’apporterait qu’un gain supplémen- 
taire de 1,4 pour 100. On peut montrer facilement 
que le gain apporté par l'emploi des harmoniques 
d'ordre de plus en plus élevé est de plus ‘en plus 
faible, quoique le rendement maximum: augmente 
en tendant vers 100 pour 100 quand on approche 
de la forme idéale (10). La signification. physique 
de l'expression (10°) est que tous les électrons sortis du 
système de commande au cours d’une période passent 
à la même époque à travers le collecteur, cette époque 
étant celle où le champ électrique du collecteur est 
retardateur et maximum en valeur absolue. 

Les auteurs reportent à un article futur la discus- 
sion des possibilités techniques de mise en appli- 
cation du principe précédent dont il est superflu 
de souligner l’important intérêt pratique pour l’aug- 
mentation de la puissance des générateurs d’ondes 
ultra-courtes. 

(À suivre.) 


ANNEXE I. 


Soit 


À 2 4 
=) = Ao+ di COS 9 + DA Sin? +... 
2Kx—0 


+ An COSAD + bhsinno +..., 


I °T 14 
Un —= = f (Rep) do 
0 & rs 
S(A—1IT » 9 


= A cos ny dÿ 


2 TR î 


5 x? | [- cup fous on nn y 
4 l Ÿ ME SA) T Ÿ 


9 SANT so 
Àj 1 1 j Sin æ 
= ——————— — ——— —n dæ |; 
T EE 2kF & æ 


k—Unr 
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de même 
SAUT 


= 

YA (UC à 
Din 2 dr 

Le S(A—A1)nT Gi 


En introduisant les fonctions 


CE , © 
. sin æ : cos à 
y == f d2, Cir = dr. 
: z j Fe 
1h vi 


a, et b, se mettent sous la forme 


2 À} 
4r2K(K— 1) 


An —= 


CRR PTETT: —1)az)—si(2knr)| 
be 2 fci(e(k—nn») — Ci(>knz)]. 


Ces coefficients tendent vers zéro quand n aug- 
mente indéfiniment comme on le voit en rempla- 
çant Ciy et siy par leurs valeurs asymptotiques 
i cosy 


Ciy = —— sy — - En utilisant les relations 
(cf. Jahnke Emde, Funktionen Tafein, D) 


À r : lé 0 
SE HYCIY— . H(JY) 


avec 


He) San 4e 
=i+ +... 


On obtient pour siy et Ciy les valeurs asymptotiques 
plus approchées 

LV PO) sim QC), 

= ; 

_ sin» P(r) + cosy Q(r) 


CIE - 
h 
avec 
FREE 
PO =IE UT RS 
À FN 
O()= = [£ — 7 +... |: 
AE 
d’où, en posant 
, 1— P(y O(1) 
f()= a AN 
un 1. REC ONES ere je NEUREE 2e Fr 
Uno LCA =nDnz) = 72 n7)| an?hÿ f' (a), 
nÀ? no 7 
bn = = [g(oCk—i)nr)—g(2knr)]  — an21j g'(x), 


« étant compris entre 2(k— 1) nr et 2knr. On peut 
prendre en particulier à = V4k(k— 1)n?r?= ny. 
En donnant à à cette valeur et en remplaçant dans 
les formules -précédentes 2, par /;, on obtient fina- 
lement 


DES 2x 5! 

En — Ji. 
nf, nt}, 
D AN 2 4! 

CE Re mec ——. 
AI, n3 À j 


Ê 


DIAMÈTRE ET VITESSE D'UN JET ISSU D'UN AJUTAGE CYLINDRIQUE 


Par Guy LITTAYE. 
Laboratoire de Mécanique des Fluides. 


Sommaire. — l’auteur mesure aux divers points d’un jet la vitesse du liquide en déterminant la 
normale en c2s points par uni procédé optique. Il observe une variation notable de la vitesse au voisi- 
nage immédiat de l’orifice, variation»conforme aux prévisions théoriques de l’auteur. 


Lorsqu'une véine liquide sort à bouche bée d’un 
ajutage cylindrique suffisamment long, on admet 
que son diamètre est égal à celui de l’orifice. Nous 
avons montré ailleurs que cette proposition ne 
pouvait être tenue pour valable dans tous les cas (1). 


droite du jet, à une distance suffisante de l’orifice, 
la vitesse est uniforme. Admettons que la vitesse 
quadratique moyenne conserve la même valeur 

la vitesse moyenne doit être plus grande pour le 
jet, l’aire d’une section droite doit être plus faible. 
Le rapport k des diamètres du jet et de l’ajutage 
doit être 0,84, ceci au moins pour une vitesse sufli- 
sante. En effet la tension superficielle : du liquide 
intervient ‘aux faibles’ vitesses pour modifier la 
valeur du diamètre de la veine liquide. La surface 
de la veine se raccorde à la surface extérieure de 
l’ajutage, sa courbure totale étant alors plus faible 
au voisinage de l’orifice qu’à une certaine distance 


(G) C. R. Acud..Se., 1940, 240, p. Goo. 


Supposons en effet la vitesse assez faible pour que 
l’écoulement soit laminaire et l’ajutage suffi- 
samment long. Dans une section droite à l’intérieur 
du tube la vitesse diminue de l’axe à la périphérie 
suivant une loi parabolique. Dans une section 


de celui-ci. Par suite de l’existence de la tension 
superficielle, ceci entraîne une augmentation de 
la pression à l’intérieur du jet, et, corrélativement, 
une diminution de la vitesse du liquide et ‘une 
augmentation du diamètre de la veine à mesure 
qu’on s'éloigne de l’orifice (?). Cet effet sera d’autant 
plus sensible que la vitesse sera plus faible. La 
valeur 0,8, se présente ainsi comme une limile 
inférieure de k, atteinte seulement aux vitesses 
élevées et pour des ajutages de longueur suffisante. 
Pour des ajutages courts, toutes les valeurs de k 


(2) A cet effet, qui est prépondérant, s'ajoute le travail 
nécessité par l’augmentation de la surface libre du liquide; 
ce travail entraîne, lui aussi, une diminution de la vitesse 
du jet. 
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comprises entre 0,84 et 1 pourront être obtenues; 
pour de faibles vitesses du jet, la valeur de X pourra 
même être largement supérieure à l'unité. Ces 
précisions ont été vérifiées par l’ensemble de nos 
mesures de diamètre des veines liquides (#). 

Nous avons cherché à vérifier directement la 
variation de vitesse du jet au sortir de l’orifice. 
Pour cela nous avons employé le dispositif suivant. 
L’ajutage, un morceau de tube capillaire horizontal 
de longueur convenable, présente une section 
elliptique dont le grand axe est vertical. La section 
droite du jet est alors une ellipse dont le grand axe 
est alternativement vertical et horizontal lorsqu'on 
s'éloigne de l’orifice. Le jet est éclairé par un faisceau 
lumineux parallèle et horizontal qui projette l'ombre 
de la veine liquide sur une surface photographique 
sensible. Deux écrans verticaux sont placés respec- 
tivement au-dessus et au-dessous du jet; eux aussi 
portent ombre sur la surface sensible. Chaque 
onde le long du jet joue le rôle d’une lentille cylindro- 
sphérique; le faisceau qui la traverse vient former 
un trait lumineux, normal au jet, sur la surface 
sensible, lorsque celle-ci se trouve à une distance 
convenable. En utilisant un orifice presque circu- 
laire, on doit placer la surface sensible à une distance 
notable. Les traits lumineux s'étendent alors, laté- 
ralement, à une grande distance de la veine liquide 
dans l’ombre portée par les écrans. Ils permettent 
de déterminer, en chaque point, le rayon de cour- 
Pure r de la trajectoire du jet et, par suite, la vitesse 
v=—\rg du liquide, g représentant l'intensité de 
la pesanteur. ; 

Les résultats obtenus par cette méthode ont, ici 
encore, vérifié nos prévisions, comme le montrent 
les exemples suivants. Le liquide utilisé est l’eau. 


I. Pour un ajutage de longueur suffisante, la 
vitesse du jet croît rapidement au sortir de l’orifice, 
ce qui correspond à une diminution du diamètre 
de la veine liquide. 


Diamètre de l’orilice, 4 = 1,00 mm; 


longueur de l’ajutage, {= 20,8 mm; PATIO": 


@) Thèse, l’aris 1942. Par exemple avec l’eau et un tube 
capillaire, de diamètre 0,84 mm et de longueur suffisante, 
les valeurs de k à diverses vitesses sont : 


Vitesse cm/s..... MEN 220 318 
LORS PR MR AE 0,93 0,90 0,87 
Avec l’eau el un orilice en paroi mince de diamètre 0,192 mm 


la valeur de k est de 1,28 à la vitesse de 1 16 em/sec. 


Ravon de courbure 
déterminé au milieu 
de chacune 
des ondes successives (cm). | 


Vitesse movenne le long 
de chacune 
des ondes successives 
(cm/s). ÿ 


Dome DEN ln DD 


Débit, 1,50 g/s; diamètre du jet, 0,89 mm;  Æ= 0,89: 

IT. Pour un ajutage de plus faible longueur. 
la répartition parabolique des vitesses, dans une 
section droite du tube, ne peut s'établir parfaitement ; 
l'effet correspondant peut être compensé par l’action 
de la tension superficielle du liquide et la vitesse 
du jet varie peu au sortir de l’orifice; ceci corres- 
pond à une variation faible du diamètre de la veine 
liquide. 


Diamètre de l’orilice, 4 = 0,98 mm: 


longueur de l’ajutage, {/ = 2,92 mm; DNA A0 
Rayon de courbure | 
déterminé au milieu 100 103  LOAMIION SE LOU 


des ondes successives (em). 


Vitesse moyenne le long 
de chacune é : RETER ; : 

: (UN DS ES DES AO PT 

des ondes successives | 


(cm/s). 


Débit, 2,29 g/s: 


diamètre du jet, 0,95 mm; Æ—=0,95. 

III. Pour de faibles vitesses, l’effet de la tension 
superficielle devient prépondérant : la vitesse du 
jet diminue rapidement au sortir de l’orifice. C’est 
ce que montre la photographie ci-contre sur laquelle 
on aperçoit nettement la diminution du rayon de 
courbure, de la première onde aux suivantes. A 
cette diminution de la vitesse du jet, correspond 


une augmentation du diamètre de la veine liquide. 


Diamètre de l’orifice, / = 0,98 mwn; 


longueur de l’ajutage, { = 2,092 mm; TES OS 


Rayon de courbure, 12,1 cm; vitesse du jet, 103 cm/s. 


Débit, 1,08 g/s; diamètre du jet, 1,12 mm; 

La méthode précédente de détermination de la 
vitesse d’un jet liquide a été appliquée avec succès 
à la mesure de la tension superficielle des liquides 
par la méthode de vibration des jets. Nous revien- 
drons ultérieurement sur les résultats obtenus. 


Manuscrit reçu le 15 mai «+943. 
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